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Abstract 

A system of a metastable phase with a pseudo continuous set of the 
heterogeneous centers is considered. An analytical theory for kinetics 
of decay of the metastable phase in such a system is constructed. The 
free energy of formation of a critical embryo is assumed to be known in 
the macroscopic approach as well as the solvatation energy. The theory 
constructed below is based on the quasistationary approximation for the 
nucleation rate. 



1 



1 Introduction 



The theory considered here will be based on the macroscopic approximation of 
the height of the activation barrier. This approximation lies in the base of the 
classical theory of nucleation. The theory presented below forms the part in the 
program announced in Q . All bibliographical remarks can be found in || . 

The further investigation of the heterogeneous condensation requires to take 
into consideration the process of the heterogeneous decay on the band of centers 
with the different activity. Speaking about the activity of the heterogeneous 
centers we mean that the different activity initiates the different height of the 
activation barrier A_F, i.e. the difference between the free energy of the critical 
embryo F c and the free energy of solvatation G. The set of the different activities 
of the heterogeneous centers is so dense that we can regard it as a pseudo 
continuous one. In the atmosphere this fact is proved experimentally. In the case 
of the " solid nucleus of condensation with the weak interaction" the continuous 
size of the nucleus ensures the continuous set of AF. 

Nevertheless the theory of the heterogeneous decay was constructed only for 
one type of the heterogeneous centers. So, the task to construct the kinetic 
theory of the decay for the system with the continuous set of the heterogeneous 
centers is rather essential. It will be completely fulfilled here. 

We shall use the ordinary physical assumptions formulated in j(| including 
the following ones: 

• the thermodynamic description of the critical embryo, 

• the random homogeneous space distribution of the heterogeneous centers, 

• the free-molecular regime of the droplets growth, 

• the homogeneous external conditions for the temperature and for the pres- 
sure, 

• rather a high activation barrier for the nucleus which are really going to be 
the centers of the droplets and aren't absolutely exhausted in the process. 

As far as the most interesting characteristics of this process are the numbers of 
the heterogeneously formed droplets of the different activities we shall estimate 
the accuracy of the theory by the error of the obtained solutions for these values. 
The whole process of the metastable decay can be split into the two periods: the 
period of the essential formation of the droplets and the period of the essential 
consumption of the metastable phase. At first we shall investigate the period 
of the the droplets formation. The unit volume is considered. 

The characteristic moment t» of the process of formation can be chosen as 
the initial moment. The values at the moment i* will be marked by the lower 
idex "*". 
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We assume the total number of the heterogeneous centers to be constant in 
time. 

We shall define the activity of the heterogeneous center as some parameter 
w which is proportional to the height of the activation barrier 

AF(w) = AF | w=0 -\w (1) 

with some positive parameter A. 

The total number of the heterogeneous centers with the given activity w will 
be marked by r) to t{w). Naturally r]tot(w) is rather a smooth function of w. 

The density of the molecules in the equilibrium vapor is marked by noo , the 
density of the molecules in the real vapor is marked by n. The power of the 
metastability will be characterized by the value of the supersaturation 



We shall define the super-critical embryos as the "droplets". Every droplet is 
described by the number of the molecules v or by the linear size p — v 1 ^ . Due 
to the free-molecular regime of vapor consumption we have 

d P , -l 
-dt =CaT 

where a is the condensation coefficient and r is some characteristic time between 
collisions obtained from the gas kinetic theory. 

The frontal type of the size spectrum allows to introduce the frontal size z 
according to 



f (ar^dt' (2) 
Jo 



Until the coalescence which isn't considered here equation (j^) ensures the 

growth of z in time and can be inverted 

*M = f™-^ (3, 

Hence, all values dependent on time become the values dependent on z and 
the relative size x = z — p can be introduced. During the whole evolution the 
droplet has one and the same value of the variable x. Consider t{x) as the 
moment when the droplet with the given x has been formed (as a droplet). We 
can see that the functions of time are the functions of x . Hence, we can see 
that the kinetic equation can be reduced to the fact that every droplet keeps 
the constant value of x. To reconstruct the picture of the evolution one must 
establish the dependencies t(z) and ((x). 
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2 Condensation equations system 



We shall mark by the argument oo the total values formed during the whole 
condensation process. 

We must take into account the reduction of £(0) to some value due to 
the consumption of the vapor molecules in the process of solvatation. For 
we have 

*. = C(0) - /dW?MwKH (4) 

where v e is the number of the molecules of the condensated substance in the 
equilibrium embryo. The discrepancy between $* and £(0) for the monotonic 
interaction between the molecules in the liquid phase and the heterogeneous 
center can be estimated as 

**-C(o) > v tot ^^ (5) 

Here vj tot is the total number of the heterogeneous centers of all types: 

v tot = J dw Vtot (w) (6) 

The following statements are valid in the further consideration: 

(1) The main role in the vapor consumption during the evolution is played 
by the super-critical embryos, i.e. by the droplets. 

(2) The quasistationary approximation for the nucleation rate is valid during 
the period of the essential formation of the droplets. 

They are proved analytically. 

The justification of the second statement uses the estimate for the times if of 
the establishing of the stationary state in the near-critical region which can be 
found in Q| (for heterogeneous barrier the consideration is quite analogous). 

Let f s be the stationary value of the size distribution for the heterogeneously 
formed droplets measured in the units of 7ioo . It can be presented in the following 
form 

fa = fdC{x),w)ri(x,w) (7) 

where r](x, w) is the number of the heterogeneous centers free from the super- 
critical embryos and is given by the following formula [El 

= iy+exp(-AF(C, W ))T 
JC TrVSA^AcKanoc 1 ' 

where W + is the number of the molecules interacting with the droplet in the 
unit of time, A e v is the width of the equilibrium distribution 

A e v= exp(—F(v) + G) 



4 



and A c v is the halfwidth of the near-critical region 

cV ~ I (dlH\ 11/2 

I V Q v 2 )V=V C | 

Index "c" marks the values for the critical embryo. Certainly, W + , A c v and 
A e v are some smooth functions of w and we neglect this dependence. 

We shall mark by noog(w) the total number of the vapor molecules in the 
heterogeneously formed droplets on the centers of the activity w. To simplify 
the formulas we shall use 



0(w) = r}{w)/r]tot{w) 
We obtain for g, 9 the following equations 

g(z,w)= (z - x) 3 f ( (((x), w)rj(x,w)dx (9) 
Jo 

9{z, w) = exp(-n OG f f ( (((x),w)dx) (10) 
Jo 

As far as we measure the accuracy of the theory in terms of the error in the 
droplets number we define these values as the following ones: 

N(z,w)=r hot (w)(l-6(z,w)) (11) 

The total number of the droplets is 

iV tot = J r] tot (w){l-e{z,w))dw = J N(z,w)dw (12) 

For the majority of the types of the heterogeneous centers the following ap- 
proximations for the nucleation rates are valid during the period of the essential 
formation of the droplets 

f ( (((x),w) = f c ($*,w) \ w=0 cxp(r (C ^l )exp( W A) (13) 

where 

r.-..^«)|^. ,14, 

and AF is the height of the heterogeneous activation barrier. The validity of 
these approximations can be analytically justified for the heterogeneous embryos 
with the interaction between the center and the molecules of the condensated 
phase weaker or equal than the function reciprocal to a space distance. Then 
we can imagine the center as the hard sphere with a weak interaction on which 
the embryo is formed. 
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The dependence of T on w is rather a weak one. So we can put 



rw = ru (15) 

for some essential part of the activity spectrum. 

Using the conservation laws for the heterogeneous centers and for the molecules 
of the condensated substance we obtain for g, 6 the following equations 

f z a tot 
g{z,w) = f* (z-x) 3 exp(-r^-)^xry tot Hexp( U ;A) = G w {g tot ,0) (16) 

dwg(z,w) (17) 
f z a tot 

6(z,w) = exp(-/» exp{\w) noc / exp(-r^— )cte) = S w {g tot ) (18) 

where /* = w = 0) 

These equations form the closed system of the equations for the condensation 
kinetics. This system will be the subject of our investigation. 

The size spectrum can be found as the following one 

f(x,w) = U exp(A W ) exp(-r I dW9 J Xl w) )9(x, w)g tot (w) (19) 



3 Iteration procedure 

The systems like (p~7|)- (p~8|) are the ordinary solved by the iteration procedure. 
It can be constructed by the following way: For the initial approximations we 
choose: 

go{z,w)= 6 = l (20) 
Ordinary the recurrent procedure is defined according to 

g i+1 (z,w)=G w (gf t ,e i (w)) (21) 
gf\z)= J dw 9l {z,w) (22) 

d i+1 (z,w) = S w (g$ ot ) (23) 

The remarkable monotonic properties of G, J dw and S ensures the guaranties 
for the error of the approximation. In this situation they can be observed also 
in the manner analogous to the case of ordinary heterogeneous condensation. 
Certainly, for the calculation of the iterations we must use some expression for 
Vtot(w)- So at first we shall consider the limit cases and obtain some estimates. 

An interesting fact goes from the absence of r\ to t in the r.h.s. of (p3|). So 
for some rather an arbitrary r) to t the power of exhaustion will be determined 



G 



only by w. The centers with the high activity w are almost exhausted during 
the process of condensation. The centers with the relatively small activity w 
remain unexhausted. The intermediate region has rather a small size. To give 
the qualitative estimates let us notice that the supersaturation £ appears only 
in the function fc. For the behavior of we can obtain the following estimate 

f c ~ /* exp(~constx t ) 3 < e < 4 (24) 

This estimate goes from the general obvious fact that the spectrum of the 
droplets must be wider than the monodisperse one and the intensity of for- 
mation must decrease in the x scale. So we see that the width Ax of the 
spectrum is a well defined value. We must stress that (^4|) is valid not only for 
the separate process of condensation on the centers with some separate activity 
but also for condensation on the centers with the set of all activities. 

The most rough approximation for the exhaustion of the heterogeneous cen- 
ters can be obtained from the first approximation in the iteration procedure. 
We have 

9 = cxp[— /» exp(Aw)n oc z] (25) 

and for the final value 

final — exp[-/* exp(Xw)n 00 Ax] (26) 
Let us define w < and iu> according to the following expressions 

w< = w - - (27) 

w> = w + ~ (28) 

where 

- 1 1 ( ln2 ) (29) 
A f^riooAx 

e - 1 

For w > w> almost all heterogeneous centers are exhausted. For w < w < 
almost all heterogeneous centers remain free. This remark allows to rewrite the 
expression for g tot in the following form 

rw> roc z 3 

g tot = dwg{z,w)+ r) tot (w)dw (30) 

J w< Jw > n oo 

The size of the intermediate region has in the scale of w of the order of 1/A. 
This size corresponds to the variation of AF of the order of 1. As far as AF 
has rather big values, the relative variation of AF in the intermediate region is 
rather small. So we can put in this region 

Vtot(w) = r?» = const (31) 
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We can spread this approximation over the region w > iu> because the real 
activity of these centers isn't important. Certainly, we must obtain the boundary 
Wmax of this region from the following equation: 

nam POO 

rj*dw = r)*(w max -«;>) = / r] tot (w)dw (32) 

For w < w < the accuracy of the approximation is not so essential because these 
centers remain free. 

After the scale transformations Tg tct /$* — > G; Tg/§* ~ ► g we have 

g = exp(Aw) ( (z - x) 3 exp(-G)0dx (33) 

G = / dioff (34) 



8 = exp[— /*rtoo exp(Aw) / exp(— G)dx] (35) 
Jo 

By the appropriate shift of w one can put w max to the zero. By the appropriate 
scale of w one can put A = 1. By the appropriate choice of z we put the 
coefficient in ( |33| ) to unity. As a result we have only one parameter - the 
coefficient in the last equation. We shall mark it by "A". 

The final form of the system of the balance equations is the following one 

g = exp(w) / (z — x) 3 exp(—G)0dx (36) 
Jo 

G = [ dwg (37) 



6 = exp[-Aexp(w) f exp{-G)dx] (38) 
Jo 

An iteration procedure can be constructed by the following way 

9i+i = exp(iu) / (z — x) 3 exp(— Gi)6idx (39) 
Jo 



Gi= / dwgi (40) 

J — oc 

9 i+ i = exp(— Aexp(w) / exp(— Gi)dx) (41) 
Jo 

go = 0o = 1 (42) 
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So, we have in the first iteration 

g 1 = exp(w)^- (43) 



z 4 



Gi = / dw 9l = — (44) 

J — oo 4 

0i = exp(-Acxp(u;)z) (45) 
The second approximation gives for 9 the following result 

f z z 4 

9 2 = cxp[— Acxp(w) / cxp(— — )dz] (46) 
Jo 4 

and for the final value 

O^aai 2 = cxp[-Aexp( U ;)1.28 1 / 4 ] (47) 

Then for N tot we have 

r° 

N tot = [1 -cxp(-AcxpH±.28 1/4 )](fo;r/» (48) 

J — oo 

For (7 in the second approximation we have 

gi = exp(ui) / (z — x) 3 cxp(— Aexp(w)x) exp(— x 4 /4)dx (49) 
Jo 

The analytical expression for these and the further iterations can not be ob- 
tained. So we shall investigate the reduced iterations. In this procedure Oi is 
substituted by in the expression for gi+\. So, we have 

6i = cxp[-,4cxp(u;)z] (50) 

gi = exp(w) / (z — x) 3 exp(— Aexp(w)x)dx (51) 
Jo 



1 -,3 ^ 2 



gi = — z exp(— w)+-^-zexp(— 2w) — cxp(— 3w)+cxp(— 3w)6A cxp(— A exp(w)z) 

(52) 

G\= I (z — x) 3 I exp(w) exp(—Aexp(w)x)dxdw (53) 

JO J-oo 

Gi = / (z-x) 3 — [1 - exp(-Az;)]dx (54) 
Jo 

After the decomposition of the last subintegral expression we have 

G 1 = z 3 I -3z 2l -±+3z I -l- I -l (55) 
AAA 
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where 

f z 1 

I Q = — [1 -exp(-Ax)]dx (56) 
Jo Ax 

h=z _l_eM z Az) (57) 

h = T + i ^exp(-Az) - i(l - exp(-Az)) (58) 

J 3 = y + ^j(^) 2 exp(-Az) + -^Az cxp(-Az) - -1(1 - exp(-As)) (59) 
For Jo we have the following result 

f Z 1 ~ 6Xp( ~^ = E 1 (z) + 7 + ln(z) (60) 

Jo y 

where 

7 = 0.57721 Ei{z) = -Ei{z) (61) 

and _Ei is the exponential integral function. We can present this result in the 
following approximate form 

l-exp(-y) ^ + 60z + 480 

y ay ^ Z z3 + 24z 2 + 180z + 480 W 

The final approximation for 9 has the following form 

6*2 = exp{— Aexp(iy) / exp(— G±)dx} (63) 
Jo 

This expression can be calculated only approximately. Notice that G increases 
faster than a; 3 . The action of the subintegral factors like the functions exp(— constx 71 ) 
approximately leads simply to the cut-off of the region of integration when n > 3 

aq>{-constx n ) w <d{consr 1/n - x) (64) 

In order to obtain the position of the end of the process of formation we must 
solve the algebraic equation 

d = 1 (65) 

We can also use the other method. From ( [54]) after the decomposition of the 
exponent one can get the series containing some non-negative powers of A. So, 
for small A one can restrict the series by several first terms and promote the 
calculations by the following method. All integrals with the combination of 
exp(— constx) and exp(— constx 2 ) can be calculated with the help of the Boyd's 
approximation 

*/2 f°__, .2w,/ tt/2 



— < exp(z 2 ) / exp(-i 2 )^ < — 

(7T - 1)2 J z ((7T 



(z 2 + 7r) 1 /2 + ( 7r _i) 2 - ^ > J z ^ ' " (( 7 r~2) 2 z 2 + 7r) 1 /2 + 2z 

(66) 
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All terms like exp(— constx n ) n = 3,4, ... can be treated as some functions like 
the Heavisaid's one. 

The value of 02 is the base for the final approximation Ntot(w). The total 
number of the droplets can be found directly from (K33) after the integration 



riot 



N tot = 77, / [1 - 9(z, w)dw] (67) 



In the second approximation we have 

f-0 



N tot _ 



where 



?7» / [1 — cxp(— _Bcxp(— w))dw] (68) 
J — 00 

B = A f exp(-Gi)dx (69) 
Jo 

is calculated earlier. Hence, for the final expression we obtain the integral which 
is already calculated in ([SO]), (|6^). 

The alternative method is the following one. At least after the remark that 
in the essential region there is exp(— Gi) < 1 we can decompose this exponent 
as a function of G\ and get the answer for J* exp(— G\)dx in the polynomial 
form. 



4 Universal solution method 

System of equations (pq) - (B8h doesn't allow the universal solution as in the 



case of the homogeneous condensation. Equation (36) hasn't in the r.s.h.. 
But the operator in the r.h.s. of ( p6|) ensures rather a rapid convergence of the 
iterations. The worst situation for the iteration procedure is when A = 0. In 
this situation we have the universal system 

g = cxp(w) / (z - x) 3 exp(-G)cfe (70) 

G= [ gdw (71) 

J — 00 

This system has the universal solution which will be marked by Go- 
This universal form of the size spectrum is shown by Figure 1. 
The condensation process is determined by the first three momentums of the 

distribution function 



lH{z) 



I x z exp{-G)8dx (72) 
Jo 



11 



and by the zero momentum 

= [ cxp(-G)6dx (73) 
Jo 

In the pseudo homogeneous case it is determined by 



fl + (z) 



f x i exp(-G)dx (74) 
Jo 

= / exp(-G)dx (75) 
Jo 



After the end of rather a short period of the intensive formation we can 
substitute in ( ]7^ ) and (|7^) oo instead of z in the region of the integration. So 
the further evolution will be determined by the first four momentums /ij(oo)Q 
As far as for G we substitute the universal solution Go (when A = 0), the values 
of Hi are the universal constants which can be obtained by the unique numerical 
solution of the last system. 

Now we return to the previous iteration procedure (with the new initial 
approximation). We can use Go as the initial approximation in the iteration 
procedure. So we have 

0i = exp[-Aexp(w) / exp(-G )dx] (76) 
Jo 

#2 = exp(u>) / (z — x) 3 exp(— Go) exp(— Aexp(w) / exp(— Go)dx')dx (77) 
Jo Jo 

/0 r-z f-x 

exp(w) / (z — x) 3 exp(— Go) exp(— A exp(w) / exp(—Go)dx')dxd'w 
-oo Jo Jo 

(78) 

The decomposition of the exponents gives 

G2 = 5(7TT)I Pl(z) (79) 

where 

(z - x) 3 exp(-G Q )r 00 (x)dx (80) 



^00(2:) 



/ exp(-G )dx (81) 
Jo 



1 The main role is played by 
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Then the expression for the function B(z = oo) will be the following one 



f°° f°° ^ (—A) 

B(z = oo)=A cxp(-G 2 )rfx = A / exp(-P (^)) exp(- V j- — T 
Jo Jo i=1 \ l + 1 



yPi(z))dz 
(82) 

The decomposition of the last exponent leads to 



i=l j=0 

where 



Cij = / e X p(-P (^))P/dx (84) 
Jo 

are the universal constants. For the total number of the droplets with the given 
activity we have 

N to t(w) = ?7*(1 - exp(-B(z = oo) exp(to))) (85) 

For N tot we repeat expressions (|67j), (|68|), ( |60[ ) with the new value of P and 
obtain the final result. 

In the same manner one can take into account the deviation of r\* from 
the constant value. After the decomposition of 77* into the Tailor's series we 
substitute the initial approximation go and act in the manner similar to ([76]) - 
(p5[). As a result we have some expansions in the powers of the parameter A, 
the derivatives of r\tot and the universal constants. 

There is another possibility to observe the universal solution in the wide 
class of the possible situations. We shall speak about the case of the "developed 
spectrum of activities" in the following situation. Earlier we extracted in the 
spectrum of activities the three regions: the region of rather an active centers 
(they are almost exhausted), the region of the centers with the small activity 
(they remain practically free during the whole process) and the intermediate 
region. We have noticed that in the majority of the situations the intermediate 
region has a relatively small size in comparison with the active region (the 
unactive region has no size - it can be spread till infinity). In this case one can 
fulfill some further simplifications. 

Let us transmit the point w = to the special activity for which 

\ w =o (00) = - 

i.e. the half of the centers became the centers of droplets. Then the system 
of the condensation equations after the rescaling of w can be written as the 
following one 

g = a exp(w) / (z — x) 3 exp(— G)9dx (86) 
Jo 
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g{w)dw (87) 

-oo 

6 = cxp(— ai cxp(w) / cxp(— G)dx) (88) 
Jo 

where w + is the upper boundary of the spectrum and an, ai are some constants. 
After the rescaling of z, x one can put a = 1. The condition of the choice of 
w = gives for oi: 

In 2 

J cxp(—G)dx 

So, there remains only one parameter w+. 

Note that the result of the process of condensation will lead to the following 
functional form for the number of the heterogeneous centers 



(90) 



/•OO 

V = Vtot cxp(— const exp(w) / exp(— G)dx) 

Jo 

after the process of decay or 

/oo 
exp(a; — G)dx) (91) 
-oo 

after condensation under the dynamic conditions where Vtot is the real initial 
(approximately constant) value of the number of the heterogeneous centers. The 
value of v* must be now substituted by Vtot- In any case as far as the integral 
gives some constant value the functional form is the following one 

V = Vtot exp(— const exp(w)) (92) 

or 

Oinit = exp(— const exp(w)) (93) 

After the observed process of condensation the cut-off of the spectrum w + 
will be unessential when 8i n i t {w + ) — > 0. This initiates the more natural defini- 
tion w + — > oo of the activity spectrum without rather an artificial parameter 
w+. 

Now we substitute this functional dependence into the system of the con- 
densation equations. Then we fulfill the integration over w and come to 



/*- 

Jo a 



(z - x) 3 exp(-G) 
+ } n2 - f" exp(-G(x'))dx r 



g - r 2 (94) 



where the coefficient a is initiated by the new form of the activity spectrum. 
For simplicity one can cancel the coefficient 

In 2 



J exp(-G{x'))dx' 



14 



As far as one can see that J Q °° exp(— G{x'))dx' has the values which are situ- 
ated not so far from 1 the cancellation of this coefficient leads only to some 
unsufficient shift of the zero point in the scale of the activities. 
Then we come to 

Jo a + J Q cxp(-G(x'))dx' 

The main essential problem is to establish the form of the size spectrum. 
Then the system of the condensation equations can be reduced to the algebraic 
system. 

Now we are going to show that the form of the spectrum has the weak 
dependence on the value of the parameter a. In the case a — > the last equation 
can be reduced to 



Jo L exp(-G(x'))dx' 



Jo 

which resembles the homogeneous case but with the increasing function J Q Z exp(—G(x'))dx' 
in the denominator. 

When a-toowe have the pure homogeneous case: 



G 



[ (z - xf cxp(-G)dx (97) 
Jo 



Recall that in the absence of this function the spectrum exp(— G) has rather 
a simple approximate qualitative behavior: 

exp(-G) - 1 z < z bound 

exp(-G) - z > z bound 

where the parameter zbound corresponds to G(zbound) = 1- 

Note that this function isn't so sharp function of z as (z — x) 3 is. The 
behavior of this function is the following 



/" 

Jo 



exp(-G(x'))dx' ~ z z < z hound 



/ exp(-G(x'))dx' - const z > z bound 
Jo 

So the qualitative behavior of G will be the same, only the value of the charac- 
teristic length of the spectrum z bound will be another. 

In order to get the same length of the spectrum we rescale z in such a manner 
that the parameter a + c appears in front of the external integral. 

Then we come to the following equation 



/ -2— A- 1 



{z - xf cxp(-G) 
• • — exp( G(x') )d.r'' 



g = a :\ v.r?>z' *» (^) 
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This equation has the following asymptotes: 



when a — > oo and 



G = [ {Z ~ Xf i?t G)dX (99) 

G=f\ {z - x)3 f eXP }-^ ) dx (100) 
Jo \ J Z exp(-G{x'))dx> 1 ' 



when a = 0. 

The value of c = 0.08 is determined in order to have the approximate coin- 
cidence of the characteristic lengths of the spectrums in these two limit cases. 
Our qualitative analysis shows that the form of the spectrum / <~ exp(— G) 
will be approximately one and the same for all values of the parameter a. The 
numerical results shown in the upper part of Figure 2 confirm this conclusion. 
We see that the curves fr a =o) an d /(o=oo) practically coincide. 

Certainly, the integral in the denominator can influence on the process in 
another manner. It can attain some value and violates the correct rescaling of 
the spectrum to the standard length. The maximum of such an influence will 
be attained at a ~ c. This case is drawn also and one can see that in this 
situation the flat region of the spectrum has the length of the flat region of the 
spectrum in the situation a — oo and the inclined region has the slope as the 
inclined region in the situation a = has. This doesn't lead as it is seen from 
this Figure to any essential deviation of the form of the spectrum. 

In the lower part of this Figure the boundaries for the spectrum from below 
fbei and from above f a b for c = 0.08 and arbitrary values of a are drawn. One 
can see that the form of the spectrum is the approximate universal function. 

One can investigate now the situation with some finite value of w+. The 
analogous transformations lead to the following equation 

G=f — — ( z x ) CX P( ^ [l-cxp(-(a+consf f exp(-G{x'))dx') exp(w+)]dx 

Jo lire + t& Jo eM-G(x'))dx> J 

(101) 

For w + = — oo one has the homogeneous situation. For w + = oo one has 
the already analyzed situation which as it has been proved lies not far from 
the homogeneous one. So, we have the same limits as the boundaries and can 
assume that here the spectrum is also near the universal one. 

Note that the physical reason of this universality is the sharp form of the 
function (z — x) 3 which induces the approximate rectangular form of the spec- 
trum. 

Now we return to the previous definitions of the parameters made before 
the integration over w. Let us explicitly extract the intermediate and the active 
regions. One can see that for the two different activities w\ and w 2 the following 
equation is valid 

= ln(0(w2,a)) (1Q2) 
exp(tyi) exp(ui 2 ) 
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The same is valid for the final values. The value of ln ( 6> ("'^)) j s invariant for 
the different activities acting in one and the same process. So, one can put 
w — w++ as the boundary between these two regions, where w++ is equal to 2 
or 3. In the same manner we can separate the region of the unactive centers by 
the boundary u>__ = — w ++ . One can neglect the substance in the droplets on 
the unactive centers and get^ 

r w + 

G = g(w)dw (103) 

J w 

We denote by n ao G+ the number of the condensing substance molecules 
in the droplets formed on the active centers. The system of the condensation 
equations can be rewritten as 

g = oq exp(w) / (z — x) 3 exp(—G)dx 



/w ++ 
g(w)dw + G + 
-OO 

9 = exp(— eti exp(ui) / exp(— G)dx) 
Jo 

Now one can write some expression for G+ 

G+ = — (w + - w ++ )z 3 (104) 
The system of the condensation equations is the following one 
g = ciq exp(w) / (z — x) 3 cxp(—G)dx 



Jo 

r w ++ „ 
G = / g(w)dw H — (w+ — w++)z 3 

J —oo n oo 

9 = exp(— a\ exp(w) / exp(— G)dx) 
Jo 

The value of is universal, but the coefficient in the term ^(w+ — w ++ )z 3 
depends on parameters. 
Defining 

G_ = / g{w)dw (105) 



2 The absence of the active and intermediate regions < w ) means that condensation 

occurs in the pseudo homogeneous way. 
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one can propose the following estimate 



G- < (w ++ — w ) 




(106) 



When 



w+ 



w ++ > w ++ - w 



one can approximately get 



G 



— {w+ - w ++ )z 3 



Then the first equation isn't necessary at all and one needn't to consider it. 
So we needn't to put ao to 1 and can use an arbitrary scale to cancel another 
constant. The system looks like 



After the rescaling of x, z one can put ^(w + — w ++ ) = 1 and the system looks 
like the universal expression: 



Note that this law takes place only in the case of the "wide spectrum" . 

The behavior of as a function of z for the different values of w are shown 
by Figure 3. 

Despite the accurate character of the methods discussed above these methods 
don't allow the clear interpretation. That's why we pay attention to the simple 
method presented below. It can be applied for rather an arbitrary conduction 
of r?* as the function of w. 

5 Monodisperse approximation 

The level of the supersaturation leading to the cut-off of the spectrum by the 
exhaustion of the substance is practically one and the same for all sorts of the 
droplets (all sort of the centers) . Let us see the droplets of what sizes play the 
leading role in this cut-off. Analyzing the subintegral expression in the equation 
for g we realize that this subintegral expression is the very sharp function of x. 
It is less than the function 



9 = cxp(— cii cxp(w) / exp(— G)dx) 
Jo 




G = -!?-(w + - w ++ )z 3 




s bel = Q(z - x)(z - xf 



(107) 
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and greater than the function 

rv w n3 / F J f*Vtot(w) e*p{*w)dw z 4 
Sab = 6(z-a:)(,z-a:) exp( — ) exp(-/* cxp(Aw)nooz) 

(108) 

As far as the interruption of the process of the droplets formation by the 
fall of the supersaturation is investigated we shall describe the pseudohomoge- 
ncous situation. Later the generalization will be presented. Let us extract the 
approximation for this function. In other words we must extract the region of 
the droplets of the sizes essential in the vapor consumption. This consumption 
in its turn is essential when 

iRiAi (109) 

where Ax is the characteristic size of the cut-off. This value can be introduced 
due to the frontal character of the back side of the spectrum. Certainly, this 
region must have the sizes rather small in comparison with Ax because all 
iteration procedures in the homogeneous decay are based on the fact that the 
droplets formed at the almost ideal supersaturation determine the process of the 
spectrum formation. As for the differential halfwidth 8\/ 2 we have the following 
expression 

5 1/2 x = (l-^)x (110) 
As for the integral halfwidth A-^/ 2 x, it is obtained from the following equation 

x 4 

N ess x 3 = /*— exp(Aw)noo (111) 

where N ess is the number of the essential droplets obtained as N ess — /* exp(\w)Ai/ 2 xn c 
which gives 

A 1/2 x=^x (112) 

and practically coincides with Si/ 2 x. So, the subintegral function s is now split 
into the essential part where x < ^ and the tail where x > ^p. We shall 
neglect the tail and use due to rather a small size of the essential region the 
monodisperse approximation for droplets formed in this region. As the result 
we obtain the approximation for g(x) 

•M-*!^ ma) 

where N{z/A) is the number of droplets appeared from x — till x = z/A. 
As far as spectrum is cut off by the supersaturation exhaustion in the frontal 
(sharp) manner the value of g is unessential before z = Ax as a small one. After 
the moment of the cut-off it is also unessential as there is no formation of the 
droplets. So instead of the previous approximation we can use 

. . N{Ax/A) 3 .^ A . 
g(z) = —i '-^z 3 (114) 
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exhaustion of the heterogeneous centers makes the subintcgral function more 
sharp and the monodisperse approximation becomes at Ax even better than 
in the pseudo homogeneous situation. But the exhaustion of the heterogeneous 
centers makes the coordinate of the supersaturation cut-off greater than Ax and 
the monodisperse approximation becomes even more better at the moment of 
the supersaturation cut-off. Certainly, we must use N(Ax/4) calculated with 
the account of the exhaustion of the heterogeneous centers. 

The concluding remarks concern the fact that we can obtain N(Ax/A) by 
the solution of the equations for the separate condensation process because we 
need the lowest cut-off length. This length is given without the cross influence 
taking into account due to the frontal character of the back side of the spectrum. 



6 Equation on the parameters 



We shall return to the initial system and introduce there a new approximation. 
The system of the condensation equations will be the following one 



, 3 T N tot (Ax/4) s 
g = cxp(w) / (z - x) J exp(- — —z f} )9dx 



G 



gdw 



r z r 

9 = exp[— Aexp(w) / cxp( 

Jo 



r N tot (Ax/4) 



z 3 )dx] 



(115) 

(116) 
(117) 



where 

N tot (Ax/A) 



N tot {Ax/A)dw N tot {Ax/A) = r) tot (l - 6{Ax/A)) (118) 



As far as we use the monodisperse approximation we have derived the behavior 
of the supersaturation and need only the parameters of such a behavior. So 
we can leave integral syste m (|115| ) - ( 118) and reformulate it in the terms of 
7Y tot (Ax/4), i.e. use only (117) - (Hi). As we have seen the value of Ax is 
defined only by the supersaturation cut-off. Hence, we have 



T N tot (Ax/A) 3 * 
ex P("7^ Z z > dz = ( 



\l/3 



•I- 



r Af* ot (Ax/4)' 



E 



where 



,1/4 

E = / exp{-z 3 )dz ks 0.24 
Jo 

We obtain the following system 

9(Ax/4) = exp[-Aexp(w){- " " 



r N tot (Ax/4) 



f^E] 



(119) 
(120) 

(121) 
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7V* ot (Ax/4H f dw{l - exp[-AexpH(^ N J { 2 x/A) ) 1/3 E])v* (122) 

The last equation is the ordinary algebraic equation with the well known meth- 
ods of the solution. The integral can be taken according to the procedure after 
©■ 

When r/tot is some function equation ( |l22j ) remains the algebraic one. The 
number of the heterogeneous centers is obtained on the base of N tot {^-) and 



the value of g is obtained due to (115) 



7 Concluding remarks 

The description of the further periods can be given with the help of the direct 
applications of the monodisperse approximation as in the case of the decay of 
the homogeneous metastable phase. We have the following equation for the 
hydrodynamic isolated element 

l^ = C = ^-^ ot — (123) 
a at noo 

where N^ ot is the total number of droplets. The last equation can be easily 
integrated as far as the r.h.s. doesn't depend on time. So there are no problems 
in the process description until the coalescence. 
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Figure 1 

Form of the spectrum of the universal solution. 
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Figure 2 

Approximate universiality of the spectrum. 
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Figure 3 

Universal dependence in the case of the "wide spectrum" . 
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